The Dynamical Mass and Three-Dimensional Orbit of HR7672B: 
A Benchmark Brown Dwarf with High Eccentricity 
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ABSTRACT 

The companion to the GOV star HR7672 directly imaged by iLiu et al. I (I2OO2I) has moved 
measurably along its orbit since the discovery epoch, making it possible to determine its dynam- 
ical properties. Originally targeted with adaptive optics because it showed a long-term radial 
velocity acceleration (trend), we have monitored this star with precise Doppler measurements 
and have now established a 24 year time baseline. The radial velocity variations show significant 
curvature (change in the acceleration) including an inflection point. We have also obtained a 
recent image of HR7672B with NIRC2 at Keck. The astrometry also shows curvature. In this 
paper, we use jointly-fitted Doppler and astrometric models to calculate the three-dimensional 
orbit and dynamical mass of the companion. The mass of the host star is determined using 
a direct radius measurement from CHARA interferometry in combination with high resolution 
spectroscopic modeling. We find that HR7672B has a highly eccentric, e = O.SOIq q^, near edge- 
i = 97.3lo;5 deg, orbit with semimajor axis, a — 18. 3"^°'^ 



on, 
is m 



Q 5 AU. The mass of the companion 



&S.7_^ iMj at the 68.2% confidence level. HR7672B thus resides near the substellar 



boundary, just below the hydrogen-fusing limit. These measurements of the companion mass 
are independent of its brightness and spectrum and establish HR7672B as a rare and precious 
"benchmark" brown dwarf with a well-determined mass, age, and metallicity essential for testing 
theoretical evolutionary models and synthetic spectral models. It is presently the only directly 
imaged L,T,Y-dwarf known to produce an RV trend around a solar-type star. 
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1. INTRODUCTION 

Brown dwarfs have complicated atmospheres. 
Unlike stars, for which it is possible to infer bulk 
physical properties from spectra alone, the emer- 
gent radiation from substell ar objects is currentl y 



not well understood (e.g.. Gushing et al.l [2OO8I) 



More than a decade following the direct detection 



partment of Physics and Astronomy, Georgia State Uni- 
versity, Atlanta, GA 30302 
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of the first L and T dwarfs dBecklin and Zuckerman nation with Doppler rad i al vel o city fRV) ob- 



19881: iNakaiima et al.' 'iggsF'ODDenheimer et al 



19951: iReb olo ct al. 1995; Basri ct al. 1996; Kirk patrick et 
1999bir and now faced with the recent detection of 
the fi r st Y dwarfs (ICushing et al. 2011 : Liu et al. 



2011 : Luhman et al.l 2011), theoretical spectral 
models are still undergoing major developments 
as they presently do not capture all of the rel- 
evant physics i nvolved in shapin g the spectra 
of co l d bodies (lAllard et al.l Il997l: iMarley et al ' 



20021: iBaraffe et al.l 120031: iBurrows et al.l l2006t 
Saumon and Marlevl[2008[ ) . 



Factors that complicate the interpretation of 
brown dwarf spectra include the necessity to si- 
multaneously model: the opacity of molecular 
species having millions of absorption lines, which 
results in an ill-defined continuum; the formation 
and settling of dust grains; non-equilibrium chem- 
istry resulting from convective mixing; and tem- 
poral changes from we ather pattern s, among other 



phen omena ( Gushi ng et al 



20081: iMarlev et al 



20101 ). Furthermore, the basic model input param- 
eters, such as mass, radius, age, metallicity, and 
effective temperature, are often degenerate with 
one another, particularly for the faintest objects 
for which only broadband photomet ry or low res- 
olution spectroscopy are a vailable ( Dupuv et a.1.1 



2009a; Janson et al. 2011: Galicher et al. 



l2Qll|). 

In order to improve our understanding of low tem- 
perature atmospheres and guide the development 
of more sophisticated models, it is necessary to 
measure one or more of these physical properties 
independently of spectra and photometry. 

Substellar companions found orbiting nearby 
stars serve as useful laboratories for calibrating 
theoretical models. Properties of the host star are 
more readily measured and may be used to in- 
fer properties of the companion, such as chemical 
composition or possibly age, under the assump- 
tion that the star and brown dwarf formed from 



the s ame material at the same time (|Pinfield et al 
20061 ). A number of these "benchmark" systems 



have recently been discovered (Lane et al. 2001: 


Close et al.l 2005: 


An et al. 


2007; Ireland et al. 


20081: IDuouv et al. 


l2009al: 


Bowler et al. 20091 


Biller et al.l2010t Wahhai et al. 2011). 



seryations ( Stassun et al.l 120061 : iJohnson et al 



It is also possible to measure the mass of 
brown dwarfs (and extrasolar planets) using or- 
bital dynamics. This has been accomplished in 
the past using the transit technique in combi- 



alibllaj ). However, substellar companions with 
large semimajor axes and proximate distances 
may be imaged directly (spatially resolved from 
their host) using adaptive opt ics (AO) a nd high 



contrast imaging technology ( Marois et al.l 1200 



Oppenheimer and Hinklevll200 



i 



ml 



Absil and Mawed 



2OIOI : ICrepp et al.l \201)L The companion sky- 
projected orbit can then be traced by measur- 
ing its position relative to the star over multiple 
epochs. This type of astrometry has been used for 
decades to ch aracterize the orbits of binary stars 
(|Aitkenlll919h and is significantly easier to realize 
in practice compared to monitoring the "wobble" 
of a star in the sky relative to an inertial reference 
grid with < 0.1 mas accuracy. 

If Doppler measurements of the primary star 
are also obtained over a comparable time baseline, 
one can construct a three-dimensional orbit. Since 
the astrometry breaks the sin(i) inclination de- 
generacy resulting from RV measurements alone, 
a Keplerian model that self-consistently combines 
both data sets provi des an estimate of t he compan- 
ion true mass (e.g.. iBoden et aLll2006l ). Knowing 
the mass of non-hydrogen-fusing objects is crucial 
becau se it governs their luminosity evolu tion in 
time (|Stevensonlli99ltlBurrows et al.lll997l ). 

It is commonly thought that the long orbital 
periods of wide-separation companions prohibits 
the calculation of dynamical masses. However, 
those with semimajor axes as large as a w 30 
AU orbiting nearby (d < 100 pc) stars may be 
characterized in several years to tens of years, be- 
cause the size of the astrometric uncertainties are 
~l-2 orders of magnitude smaller than the size 
of the orbit on the sky (thanks in part to the 
recent miniaturization in the physical size of de- 
tector pixels). Orbital solutions converge when 
both the astrometry and RVs become "unique", 
which generally corresponds to showing curvature 
or a change in the acceleration (also known as 
the "jerk")Q With accuracies of ~ 10 mas (e.g.. 



^ See ILu et aO l|200g| ) for analysis of an equivalent problem 
involving stars orbiting the supermassive black hole in the 
center of the Galaxy. The orbits are well-characterized even 
though the periods are of order thousands of years. Like- 
wise, the orbits of comets residing in the outer solar system 
may be determined with a time b aseline of only several days 
l lBernstein and Khushalanill200Cll V 



2 



the size of a pixel in the NIRC2 camera at Keck 
- PI: Keith Matthews), it is possible to calcu- 
late the inclination, as well as other parameters, 
in only several tenths of an orbital phase wrap 
(tens of degrees change in the true anomaly) de- 
pendin g on the orbit qrienta t ion and observing 
dates ( Konopackv et al. 2010l : Currie et al.l 12011 : 
iDupuy and Liul l201l[ ). The constraints become 
tighter with time as the companion executes its 
motion along the ellipse. 

A natural method to identify promising high- 
contrast imaging targets is to select stars that 
show long-term Doppler accelerations (trends). 
This approach is convenient because RV programs 
have now established time baselines exceeding a 
decade. For example, the L4.5±1.5 dwarf compan- 
ion to HR 7672 (=G1 779, =HD190406) was dis- 
covered by iLiu et al using this technique. 
This GOV star was originally chosen for AO obser- 
vations because it showed unambiguous evidence 
for the existence of a distant companion with sub- 
stellar minimum mass: an « -24 m/s/yr accelera- 
tion over a 14 year baseline. HR7672B is currently 
the only directly imaged L,T,Y-dwarf known to 
produce an RV trend around a solar-type star. 

In this paper, we present an updated RV time 
series that shows significant orbital motion of the 
companion since the epoch of the direct imaging 
discovery. We have also obtained a new astromet- 
ric measurement with NIRC2 at Keck. Combining 
our observations with astrometry from the liter- 
ature, and performing a joint-fit analysis to the 
full data set, we construct a three-dimensional or- 
bit and calculate the dynamical mass of IIR7672B 
with a fractional error of only 4% using Newto- 
nian dynamics (95% confidence level). By ex- 
plicitly connecting the mass of a brown dwarf to 
its radiative properties, the results may be used 
to inform both theoretical evolutionary models 
and synthetic spectral models. This study is af- 
forded by an equally detailed characterization of 
HR7672A, for which we have modeled its high res- 
olution spectrum to determine its surface gravity 
and metallicity, acquired a spectral energy dis- 
tribution (SED) to determine its luminosity, and 
measured its radius directly using interferometry. 



2. OBSERVATIONS 

2.1. Doppler Measurements 

HR7672 was initially observed with the 0.6m 
Coude Auxiliary Telescope (CAT) at Lick Ob- 
servatory on September 9, 1987 UT. Several 
years of observations using the CAT as well as 
the Lick 3m (Hamilton) - both of which feed 
the Coude echelle spectrograph - revealed the 



star to have a linear RV trend (|Cumming et al 
1999h . We also began Doppler monitoring of 



this star at Keck on June 2, 1997 UT using the 
Hlgh-Resolutioii Echelle Spectrometer (HIRES; 



Vogt et al.l 119941 ). For each instrument, we used 



the iodine cell referencing method to measure pre- 



cise RVs (jMarcv and Butleri 119921 : iButler et al 
19961 ). Both data sets sho wed a con si stent long- 
term acceleration. When iLiu et al.l ([2002) tar- 
geted this star with AO, finding that indeed a 
low-mass companion with a wide orbit was re- 
sponsible for the trend, there was very little, if 
any, perceptible change in the acceleration. 

We have continued to monitor HR7672 at both 
Lick and Keck Observatories and have now estab- 
lished a 24 year time baseline. Our most recent RV 
data reveal significant curvature, making it possi- 
ble to place tight constraints on the companion 
orbit. When combined with astrometry, this in- 
formation is sufficient to converge to a unique or- 
bital solution and companion mass, despite having 
coverage over a fraction of a single orbit cycle, be- 
cause the signal-to-noise ratio in both data sets is 
high. 

2.2. Astrometry from High-Contrast Imag- 
ing 

HR7672B has been imaged se veral times since 
i ts di scovery. In July 2002, iBoccaletti et al. 
( 2003) detected the companion with PHARO 
( Havward et al.l 2001 ) at Palomar in each of the 
J, H , Kg bands to obtain color inform ation and 
validate the spectral-type assigned by iLiu et al 



(|2002D from K-band s pectro scopy. In September 



2006, ISerabvn et al 



( 20091) observed HR7672B 
at Palomar to demonstrate the feasibility of new 
coronagraphic technologies and the use of "ex- 
treme" AO. Both studies reported astrometric 
measurements of the position of HR7672B rela- 
tive to the primary. We have also retrieved VLT 
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Fig. 1.— Direct image of HR7672B in the K' filter taken on May 15, 2011 UT with NIRC2 at Keck, (left) 
Single image prior to PSF subtraction. The companion is located inside of the quasi-static speckle halo though 
discernible in individual frames, (right) Processed image using the ADI technique and LOCI algorithm. The 
companion is detected at 36(t with 450 seconds of integration time. The inset subimage displays a cutout 
of the central region when our speckle suppression algorithm uses all available ADI frames but avoids PSF 
subtraction over the coronagraphic spot. The stellar Airy pattern is clearly visible facilitating astrometric 
measurements. We use data sets both with and without PSF subtraction to calculate the companion position 
relative to the star. This most recent observation increases the astrometric time baseline by 3.6 years. 



archival images of the system from September 
2007 (PI: A. Boccaletti - program ID 279. C- 
5052(A)). These observations show clear orbital 
motion in a clockwise directio n (north u p , east 



(|2002l) 



left) when combined with the iLiu et al 
results. 

We obtained additional images of HR7672B on 
May 15, 2011 using the Keck AO system with 
NIRC2. Given the large flux ratio between the 
primary and companion, we used the 300 mas di- 
ameter coronagraphic occulting spot to permit rel- 
atively long exposures without saturating the de- 
tector. The spot is partially transmissive and thus 
allows us to accurately measure the companion 
separation and position angle. Our initial data set 
consisted of 20 frames recorded using the narrow 
camera in position- angle mode (15 with the K' fil- 
ter and 5 with the H filter). Although difficult to 
identify in a single exposure without prior knowl- 
edge of its approximate location, the majority of 
our raw frames did reveal the companion (Fig. 1). 

To further isolate the signal of IIR7672 B 



from residual scattered starlight, we also ob- 
tained a 30 frame sequence with the K' filter 
using t he angular different ial imaging (ADI) tech- 
nique ( Marois et al. 20061) . The parallactic angle 
changed by 12.7° during this second sequence of 
images, providing sufficient angular diversity to 
remove the spec kles and keep l ight from the com- 



panion (L afreniere et al.l 120071 ). Fig. 1 displays 



our high-contrast images before and after speckle 
suppression. 

We measured the position of HR7672B relative 
to the primary star in each data set {K' with- 
out ADI, H without ADI, and K' with ADI). 
The companion and stellar PSF core were fitted 
with a Gaussian function using least-squares it- 
eration to find the astrometric centroid position 
in each processed frame. We correct for differen- 
tial geometric distortion (warping) using the pub- 
licly available code provided by the Keck NIRC2 
astrometry support webpage. Adopting a plate 
scale value of 9.963 ± 0.006 mas pix^^ and instru- 
ment orientation relative to the sky of 0.13°±0.02° 
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Year 


JD-2,440,000 


p (mas) 


PA (°) 


Instrument 


Reference 


2001.64 


12144 


786 ±6 


157.9 ±0.5 


NIRC2 


Liu et al. (2002) 


2001.94 


12254 


794 ±5 


157.3 ±0.6 


NIRC2 


Liu et al. C2002) 


2002.54 


12473 


788 ±6 


156.6 ±0.9 


PHARO 


Boccaletti et al. C2003) 


2006.69 


13989 


750 ± 80 


155.0 ±5.0 


PHARO 


Serabvn et al. (2009) 


2007.73 


14367 


742 ± 35 


151. 8± 2.9 


NACO 


Program 279.C-5052(A) 


2011.37 


15697 


519 ±6 


147.1 ±0.5 


NIRC2 


present study 



Table 1: Astrometry measurements used for orbital analysis. 



east of north, as measured bv lGhez et al.l (|2008l ). 
we find a companion separation and position an- 
gle of 519 ± 6 mas and 147.1 ± 0.5° respectively!! 
HR7672B orbits in a clock-wise direction and is 
now significantly closer to its host star compared 
to the discovery epoch. 

To estimate the uncertainty in our measure- 
ments, we calculate the standard deviation in the 
position resulting from each data set. We include 
a characteristic 5 mas error to account for astro- 
metric bias introduced by the coronagraphic spot 
(Konopacky 2011, private communication). This 
error is added in quadrature with the uncertainty 
from the plate scale and instrument orientation. 
The results were combined using a weighted aver- 
age to determine the final uncertainties. We find 
that the measurements from each data set are con- 
sistent with one another to within the {la) error 
bars. Simulated companions were injected into the 
images and their positions recovered to validate 
our results. Table [T] lists the separation and po- 
sition aiigles used for our three-dimensional orbit 
analysis |3 

2.3. Physical Properties of the Host Star 

The mass of the companion is tied to the 
mass of the primary star through the RV semi- 
amplitude and Kepler's equation. Astrometry de- 
termines the system total mass, and Doppler mea- 
surements can break the degeneracy between the 
system total mass and individual masses. How- 
ever, we only have precise Doppler information 



^This result is consistent with the values (9.94 ± 0.05 mas 
pix~^ and 0.0 ±0.5°) provided by the observatory. NIRC2 
was not opened nor temperature cycled between 2003-2011, 
including the dates of our imaging observations (R. Camp- 
bell 2011, private communication). 

''We do not include the original Gemini imaging data from 
June 2001, which suffers from syste matic errors and a low 
signal-to- noise ratio as described in I Liu et al ] ll2002h . 



for HR7672A, since the companion is significantly 
fainter than the primary. It is therefore important 
to reliably determine the mass of the host star 
independent of dynamical considerations. 

Fortunately, HR7672A is a nearby {d = 17.77± 
0.11 pc) solar- type (GOV) star, making it possi- 
ble to: (i) measure its radius directly using inter- 
ferometry; and (ii) determine its surface gravity 
from high resolution spectra using models that are 
well-calibrated from observations of the Sun. We 
have spatially resolved the surface of HR7672 A us- 
ing the Georgia State University Center for High 
Angular Resolution astronom y Array (CHARA) 



( ten Brummelaar et al. 20051 ). We have also ob- 



tained non-iodine (template) spectra of HR7672A 
for which to model. 

In the following, we calculate the mass of 
HR7672A using an iterative version of the Spec- 
troscopy Made Easy (SME) spec tral modeling 
routine (IValenti and Fischer 2005 ) , along wit h 



Yonsei-Yale isochrones (|Demaraue et al.l |200J) 



The models self-consistently incorporate the high- 
resolution stellar spectrum measured with HIRES, 
direct stellar radius measured with CHARA, and 
stellar luminosity determined from a SED. Re- 
sults for the physical properties of HR7672A, in- 
cluding its radius, luminosity, effective temper- 
ature, metallicity, mass, and age are shown in 
Table [2j Our mass estimate is confirmed using a 
separate analysis that is model independent and 
relies only upon the empirica l mass-radius rela- 
tions from [Torres et al.l (|2010l) . 



2.3.1. Stellar Diameter from Interferometry 

HR7672 was observed on the nights of Septem- 
ber 18, 19, 2010, August 21, 2011, and October 
2,3, 2011 with CHARA at Mount Wilson. We 
used the CHARA classic beam combiner in H- 
band {Xcentrai = 1.67 /im) with the longest avail- 
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Fig. 2. — Calibrated visibility measurements and 
limb-darkened angular diameter fit for HR7672A. 
Dotted lines indicate the 3cr uncertainty level in 
the overall fit to the Bessell function. The bottom 
panel shows the residuals (Xr = 0.5). Interfero- 
metric measurements of the radius of HR7672A 
allow us to place tight constraints on its mass. 

able baselines, SlEl (maximum 330 m) and ElWl 
(maximum 313 m). Projected baseline lengths 
ranged between 245 m and 325 m during our ob- 
servations. 

Our o bserving strategy i s analo gous to that em- 
ployed in lvon Braun et al.l (|2011al) and lvon Braun et al 
(j2011b( ). We briefly repeat the principle compo- 
nents below. Our interferometric measurements 
include the common technique of taking brack- 
eted sequences of the object with calibrator stars 
to characterize and eliminate effects from tempo- 
ral variations of the atmosphere, telescope, and 
instrument upon our calculation of interferomet- 
ric visibilities. We alternate between two nearby 
calibrators, HD187923 and HD192425, during the 
observations to minimize systematic errors. 

The uniform disk, ^uQj a-nd limb-darkening 
corrected, 6'ld, diamctcrto are found by fitting 
our calibrated visibility measurements to the 
respective functions for each relation. Specifi- 
cally, we use a linear combination of Bessel func- 



tions ( Hanburv Brown et al.lll974i iBovaiian et al 



Fig. 3. — Spectral energy distribution fit for 
HR7672 (Xr = 0-17)- The (blue) sp ectrum is a 
GOV spectral template ()Pickleslll998l ). The (red) 
crosses indicate photometry values from the lit- 
erature. "Error bars" in the horizontal direction 
represent bandwidths of the filters used. The 
(black) X-shaped symbols show the fiux value of 
the spectral template integrated over the filter 
transmission. The lower panel displays the resid- 
uals around the fit in fractional fiux units of pho- 
tometric uncertainty. 



Limb dar kening coefficients were taken 



fromlClaretl (|2000[ ). The data and fit for 6'ld are 
shown in Fig. [31 The key to understanding the 
uniqueness of the solution is to recognize that 
the visibility must approach unity as the base- 
line approaches zero, and that degeneracies in the 
slope are broken by the absolute visibility value 
(as opposed to relative value) for a given base- 
line, so long as the location along the visibility 
curve is not mistaken for side-lobes in the Bessel 
function (which effectively corresponds to knowing 
the distance and luminosity class of the star) . Our 
interferometric measurements yield uniform disk 
and limb-darkening corrected angular diameters of 
6»uD = 0.567 ±0.010 mas and 6Ild = 0.584 ±0.010 
mas respectively. Combined with the Hipparcos 
trigo nometric parallax v alue of tt = 56.28 ± 0.35 



mas 



( van LeeuwenI [2007 ) ■ we calculate the linear 
radius of HR7672A to be i? = 1.115 ± 0.0217?© 
(Tabled]). 



*The limb-darkening corrected Sld corresponds to the an- 
gular diameter of the Rosseland, or mean, radiating surface 
of the star. 
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Fig. 4. — Comparison between a HIRES spectrum 
(black) and our best-fit SME model (maroon). We 
use these results along with Yonsei-Yale isochrones 
to derive the physical properties of HR7672A, self- 
consistently taking into account independent mea- 
surements of the stellar radius from interferometry 
and luminosity from an SED. 



2.3.2. 



Stellar Luminosity and Effective Temper- 
ature 



Foh owing the procedure outlined in Ivan Belle et al.l 
(|2007l ). we produce a fit to HR7672 ' s SEP based 
on the spectral templates of Pickled ( 1998f) to lit- 
erature photometry published in the references 
shown in Table [T) A GOV spectral-type provides 
the best fitting template. Interstellar extinction is 
a free parameter in the fitting process and is calcu- 
lated to be Av = 0.077 ± 0.017 mag. The SED fit 
for IIR7672 along with its residuals are shown in 
Fig. [3l We calculate HR7672's stellar bolometric 
flux to be Fbol = (1.360±0.028) x lO"'^ erg cm"^ 
s"^, and its luminosity to be L = 1.338±0. 032^0. 
Using a rewritten version of the Stefan-Boltzmann 
Law, 

Teff(K) = 2341(FbolKd)^ (1) 

where Fbol is in units of 10^^ erg cm^^ s^^ and 
0LD is in units of mas, we find IIR7672 has an 
effective temperature of Tgg- — 5883 ± 59 K. 

2.3.3. Stellar Surface Gravity, Metallicity, Mass 
and Age 

Informed by the above results, we use high res- 
olution spectroscopy, spectral modeling, and stel- 



lar isochrones to determine the remaining physi- 
cal properties of the host star. A template spec- 
trum of HR7672 was obtained with Keck/HIRES 
(A/AA = 55, 000) on Sept. 8, 2008. The iodine gas 
cell was removed from the optical path to isolate 
the absorption lines of the star. With an integra- 
tion time of 17 seconds, we achieved a signal to 
noise ratio of 280 per pixel at 550 nm. 

We analyzed the stellar spectrum using the 
LTE spectral synthesis code Spectroscopy Made 
Easy (SME) described in IValenti and Fischer 
( 20051 ). Recent versions of SME include an it- 
erative scheme that ensures consistency between 
the surface gravity derived fr om synthetic spec - 
tra and evolutionary models f|Valenti et al. 2009f) . 
The stellar spectrum is fit using an initial guess for 
the surface gravity while the effective temperature, 
metallicity, macroscopic turbulence, and projected 
rotational velocity are varied to minimize residu- 
als. With each iteration, result s are corapare d 
to Yonsei-Yale isochrones tPemargue et al.ll2004l) . 
Model grids are interpolated to yield an estimate 
for the stellar mass, radius, ag e, and composition 



(see Fig. 1 of IValenti et al.l j2009l ). The spec 



tral fitting procedure is repeated using the new 
(isochronal) surface gravity estimate in the fol- 
lowing iteration until the procedure converges - 
i.e., until the surface gravities agree. 

In the case of HR 7672, we fit the stellar spec- 
trum by holding T^ff fixed, setting it equal to 
discrete values within the allowable range derived 
from the Stefan-Boltzmann relation (see ^12.3.21) . 
Specifically, we analyzed the results for three sep- 
arate cases: T^ff = Tq, T^ff = Tq + AT, and 
Teff = To - AT, where Tq = 5883 K and AT = 59 
K f ii2.3.2| . The remaining parameters were varied 
as usualO We find that the T^/ / = Tq -f AT case 
provides the best fit to the HIRES spectrum. It 
is also the only solution for which the luminosity 
and radius match the directly measured values to 
within Icr (the other cases are consistent to within 
2a). 

With this technique, we find a stellar surface 
gravity, log g — 4.42 ± 0.06cm/s^, metallicity. 



^Alternatively, we could set the luminosity and radius input 
to SME equal to our measured values. By allowing the lu- 
minosity and radius to vary, we retain covariance between 
the parameters and demonstrate that the models naturally 
reproduce measured values while also accounting for uncer- 
tainty in the measurements. 
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[Fe/H] = 0.05 ± 0.07, mass, = 1.08 ± O.O4M0, 
and age, tage = 2.5 ± 1.8 Gyr. Fig. H] shows a 
comparison between the HIRES spectrum and our 
best synthetic model. HR7672's physical proper- 
ties are listed in (Tab. [5]). Uncertainty in the stel- 
lar mass includes propagation of errors from the 
uncertainty in Tf,ff. We adopt this stellar mass 
for all dynamical calculations. 

To further validate our estimate of the host 
star mass, we use a separate technique that is 
independent of theoretical s pectra l and evolution- 
ary models. Torres et al. ( 2010[ ) have compiled 



results for dozens of eclipsing binary stars and 
have constructed a polynomial that relates their 
measured masses and radii empirically. Using 
this relation along with the other parameters of 
HR7672A listed in Tabled we find that HR7672A 
has a mass of Af, = 1.10± 0.04Mq. This resuh is 
in agreement with our analysis based on SME and 



Yonsei-Yale isochrones. Further, the [Torres et al 



(I2OIOI) relations contain a known bias that system- 
atically over-estimates near-solar [M^, k, I.OMq) 
stellar masses by as much as « 5%. Taking this 
into account indicates that our above result of 
= 1.08 ± O.O4M0 is accurate. Finally, we 
note that our age calc ulation a ffirms the original 
1-3 Gyr estimate from iLiu et al . (2002), which is 
based on an independent analysis using six differ- 
ent diagnostics, including: comparison of absolute 
visual magnitude to the location of the main- 
sequence, rotation period. X-ray emission, Ca II 
H-l-K emission, kinematics, and lithium absorp- 
tion. 



3. JOINT-FIT ORBITAL ANALYSIS 

We use Bayesian inference and Markov Chain 
Monte Carlo (MCMC) simulations to calculate the 
companion true mass, all six orbital elements, and 
their uncertainties. The RV and astrometry data 
are simultaneously fit with a Keplerian orbit. Stel- 
lar RV variations, Vm{ti)^ are modeled according 
to: 



Kn(^j) — K [cos(cl; -I- f{tj)) + e cos(a;)] , 
where. 



K = 



'2t:G 



1/3 



m sin(i) 



\/T~~i2(M, +^)2/3' 



(2) 



(3) 



is the RV semi-amplitude, uj is the argument of 
periastron, f{tj) is the true anomaly at epoch <j, 
e is the eccentricity, i is the inclination, P is the 
orbital period, is the mass of the star, and m is 
the mass of the companion. Astrometry operates 
in the orthogonal direction and relates the semi- 
major axis, a, to the orbital period and system 
mass through Kepler's equation. 



p2 



2^3 



(4) 



such that the orbit corresponds to the sky- 
projected separation of the companion at a given 
epoch (see Equ. [5]) and system parallax. 

The following physical parameters are used as 
variables in the analysis: the companion period 
(P), eccentricity (e), inclination (i), argument of 
periastron (w), longitude of the ascending node 
(il), time of periastron passage (ip), and true mass 
(m). Several nuisance parameters are also re- 
quired, including differential RV offsets between 
HIRES and the Hamilton echelle at Lick, to ac- 
count for the different instrument settings, as well 
as RV "jitter" to account for astrophysical noise 
and ensure that the data are weighted properly 
( Isaacson and Fischer 20101 ). RV jitter was added 
directly to the measurement uncertainties, which 
results in a reduced chi-squared statistic, v'xf, 
near unity. Given the long time baseline of the 
observations, a total of six different instrument off- 
sets and seven different jitter terms were required 
(five from Lick and two from Keck; see Tab.[S]). As 
expected, the measurement uncertainties decrease 
with time as well as the relative size of the re- 
quired jitter, meaning that the RV precision and 
our ability to estimate the RV precision has im- 
proved significantly since September 1987. 

We use the Metropolis-Hastings algorithm to ef- 
ficiently explore the multi-dimensional parameter 
space. A likelihood function, L, self-consistently 
relates the data, model parameters, and prior in- 
formation, toJhe^gosteriOT^grobaH^ dis- 



tribution (|Fordll2006l: Ijohnson et"alll2011bt ). The 



likelihood function used for calculations is given 
by: 



Nrv 



AV{t, 



o-j + Si _ 



Nas 
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Fig. 5. — RV data and best-fit model. The Doppler signal shows significant curvature indicative of a high 
eccentricity, (left) RV time series showing a full predicted orbit cycle, (right) A closer view of the 24 year 
data set. 
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observations compared to the other studies. The best-fit model passes through the (Icr) errorbar of 
each measurement. 
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, NAst 

4e 



+ 



(5) 



where AV{tj) — Vj ~ Vm{tj) represents the dif- 
ference between the RV data and the Keplerian 
model RV at epoch tj; AX{tk) and AY{tk) repre- 
sent the difference in the east and north position of 
the companion between the astrometric data and 
the Keplerian model astrometry at epoch tk re- 
spectively; ai and , o'Vk are the individual RV 
and astrometric uncertainties; and si is the RV jit- 
ter for instrument setting £ (summation over the 
various instrument settings is implicit). We as- 
sume that both the RV and astrometric measure- 
ments follow Gaussian distributions. 

We randomly select one parameter per MCMC 
iteration to alter. Candidate transition proba- 
bility functions follow a normal distribution with 
adjustable width and are centered on the most 
recently accepted parameter value. We do not 
change the transition function widths following the 
"burn-in" stage. The algorithm accounts for the 
covariance between w and tp by taking steps in 
wib/o, where /o is the true anomaly of the compan- 
ion at the first RV observing epoch. This approach 
improves the MCM C acceptan ce rate and accel- 
erates convergence (jFordl 120061 ). We use uniform 
priors for each parameter. The range of values is, 
however, truncated to reasonable limits. For ex- 
ample, we only consider companion masses in the 
range 20 < m/Mj < 120. The eccentricity and 
other orbit parameters span the full range of pos- 
sible values. 

To identify the global minimum, we compare 
the results of multiple chains that explore the 
likelihood manifold starting from different initial 
states. Convergence is reached once the Gelman- 
Rubin statistical criterion is met (Gelman & Ru- 
bin 1992). Specifically, we require that: 



R{z) = 



IS{z) + M{z) 



S{z) 



<1.1, 



(6) 



for each parameter z, where S{z) is the mean of 
the variance of the MCMC chains, and M{z) is 
the varian ce of the mean of the MCMC chains (see 
Fordll2006l for details). Following convergence, the 
individual chains are combined (linked together) 
to create the final parameter distributions. We 
find that w 10^ iterations are required to provide 



a sufficiently dense sampling of the posterior dis- 
tribution. 

Uncertainty in the distance to HR7672 is self- 
consistently folded into the analysis by drawing 
random distance values from a normal distribution 
centered on the Hipparcos result for each MCMC 
iteration. The width of the distribution is set to 
match the measurement error. We account for 
uncertainty in the mass of the host star in the 
same way as the uncertainty in the distance, by 
drawing random stellar mass values from a nor- 
mal distribution centered on the stellar mass esti- 
mate (i i2.3.3p . The resuhs arc then combined with 
the MCMC chains in accordance with the above 
equations. This technique of incorporating uncer- 
tainty in the distance to HR7672 and the mass of 
the primary is robust because it preserves the co- 
variance between orbital parameters captured by 
the MCMC calculations. 

4. RESULTS 

The RV and astrometry data used in combi- 
nation with the distance measurement and stel- 
lar mass estimate allow the MCMC simulations 
to consistently converge for each parameter. Only 
a narrow range of orbital solutions exist that can 
satisfy the constraints imposed by both data sets 
simultaneously. Figures 5 and 6 display the final 
RV and astrometry data along with the best-fit 
model. The orbit well-matches the RV observa- 
tions and also falls within the (Icr) uncertainty of 
each astrometric measurement. We find a reduced 



chi-square, Xr, value of 



y2 _ 



0.96 per degree of 



freedom overall. Fig. 7 shows the RV residuals fol- 
lowing subtraction of our best-fit model from the 
data. 

The final MCMC distributions are shown in 
Fig. 8 and the physical properties of HR7672B 
are summarized in Table[3l We find that HR7672B 
has a high eccentricity, e = 0.50^Q gi[, and a near 
edge-on orbit, i = 97.3lo;5 deg. It reached great- 
est elongatio n (south-east po rtion of orbit) just 
following the iLiu et aL ( 2002 ) direct imaging dis- 
covery. The velocity of HR7672B is now increasing 
sharply as it approaches periastron. Our predic- 
tions for its location on Sept. 1 in the years 2020, 
2030, 2040, 2050, and 2060 are shown in Fig. 6. 
With an impact parameter of 67 mas, the com- 
panion will soon disappear behind the host star. 
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Fig. 7. — Doppler RV residuals. Precision RV 
techniques have improved dramatically over the 
past two decades. In addition to a long time base- 
line, HR7672 also has many samples per time be- 
cause it is part of the NASA-UC 77-Earth program. 
Some residual periodicities are noticeable but cur- 
rently insufficient to claim evidence for the exis- 
tence of an additional low-mass body with a com- 
paratively short orbital period. 



Fiher 


HR 7672 A 


HR 7672 B 


V 


5.80 




J 


4.69 


« 14.39 


H 


4.43 


14.04 ±0.14 


K 


4.49 


= 13.04 ±0.10 



Table 4: Apparent magnitudes of the HR7672 sys- 
tem. Values for HR 7672 A are from SIMBAD. 



Values for HR 7672 B are from Boccaletti et al 



(|2003D 



much like the extrasolar planet Beta Pictoris b 



has evaded detection in the past (jFitzgerald et al 
20091 : iLagrange et aLllioToj ). 



We find that HR7672B has a mass of to = 
68.7^3^Mj and thus resides in an interesting 
regime that borders, but lies just beneath, the 
hydrogen- fusing limit (Fig. 9). Only 9.9% of 
the distribution lies above the canonical value of 
K.72Mj often taken as the dividing line betwee n 



20001) . 



brown dwarfs and stars (jChabrier et al. 
We therefore conclude that HR7672B is a brown 
dwarf at 90.1% confidence (Fig. 11), assuming 
the companion has the same near-solar metal- 



licity as the primary, which we measure to be 
[Fe/H]= 0.05 ±0.07 dex. 

Figure 9 shows the companion mass distribu- 
tion before and after taking into account the un- 
certainty in parallax and stellar mass. Uncertainty 
in the companion mass is limited primarily by 
having only a partial orbit, motivating the need 
for continued Doppler and astrometric monitor- 
ing. Uncertainty in parallax is the second largest 
effect. Significant improvements in parallax over 
the current Hipparcos measurement will likely re- 
quire dedicated observations from space. Uncer- 
tainty in the mass of the primary star is a com- 
paratively small effect, since the companion mass 
depends weakly upon star mass with Doppler ob- 
servations. For example, an error of ±O.13M0 
in the host star mass would be required to shift 
the companion mass exterior to its current 95.4% 
confidence interval. Covariance matrices between 
the six orbit parameters and companion mass are 
shown in Fig. 10. 

The final companion mass distribution is nar- 
row compared to the range of possible masses 
resulting f rom s pectrophotometry. For example. 



Liu et al.l ()2002f) find that the companion has a 
mass between 55 — 78Mj, using direct K-band 
spectroscopy, a numbe r of age diagnos t ics, an d the 



theoretical models of iBurrows et al. (1997) and 



Chabrier et ahl (|2000l ). iBoccaletti et al.l (|2003l ) 



have obtained complementary photometry in the 
J, H bands to refine this estima te. Us i ng th e 



same age range of 1-3 Gyr from lLiu et al.l (|2002[ ). 
they find that HR7672B has a mass between 58 — 
llM.j. Performing: similar calculations with the 
iBaraffe et all l|200 3') evolutionary models, we find 
that HR7672B sho uld have a ma s s in t he 57 — 
74Mj range. The iBaraffe et all (|20Q3f ) models 



are in excellent agre e ment with the more recent 
Saumon and MarlevI (2008) models for L-dwarfs 



older than several Myr. Thus, in each case our cal- 
culations using orbital dynamics are already more 
precise and accurate than estimates based on ana- 
lyzing light received directly from the companion 

(Fig. m. 

5. SUMMARY & DISCUSSION 

We present the first three-dimensional orbit 
solution and mass determination for a directly 
imaged L-dwarf companion to a solar-type star. 
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HR7672A Value Technique 

0LD (mas) 0.584 ±0.010 interferometry 

Radius (i?0) .... 1.115 ±0.021 interferometry, parallax 

Luminosity (^0) 1.338 ± 0.032 spectral energy distribution 

Tgff (K) 5883 ± 59 Stefan-Boltzmann relation 

log[i?^^] -4.854 ± 0.025 HIRES spectroscopy 

^macro (km/s) . . 4.8 iterative SME 

ysini (km/s) .. 2.1 ±0.7 iterative SME 

[Fe/H] 0.05 ±0.07 iterative SME 

log g (cm/s^) . . . 4.42 ± 0.06 iterative SME 

Mass (Mq) 1.08 ±0.04 iterative SME 

Age (Gyr) 2.5 ±1.8 iterative SME 



Table 2: The physical properties of HR7672A. We measure the stellar radius directly using CHARA in- 
terferometry. The luminosity is found by fitting a spectral energy distribution from photometry obtained 
over a broad wavelength range. The effective temperature derived from the Stefan-Boltzmann equation is 
held fixed as input to SME (sampled separately at T^f / — 5883 ±59 K) . Several SME iterations result in 
convergence of the stellar surface gravity (log g) between the spectral model and Yonsei-Yale isochrones. 
The final estimated mass, age, and radius are consistent with the measured stellar luminosity and radius 
to within la . Our derived stellar mass is also consistent with the model-independent empirical r elations of 
Torres et al. I (I2OIOI) . Our derived stellar age is consistent with the six different age diagnostics of iLiu et al 



(|2002n . HR7672A is a GOV star and thus a convenient calibrator. Its physical properties may be used to 
infer those of IIR7672B, such as metal content and age. 



HR7672B 


weighted mean 


68.2% CI 


95.4% CI 


mass (Mj) 


68.7 


65.6-71.1 


63.2-74.6 


P (yr) .... 


73.3 


70.4-75.5 


68.5-79.0 


a (AU) . . . 


18.3 


17.8-18.7 


17.5-19.3 


e 


0.50 


0.49-0.51 


0.48-0.52 


i (deg) .... 


97.3 


96.8-97.7 


96.4-98.1 


w (deg) . . . 


259 


257-261 


254-263 


n (deg) . . . 


61.0 


60.6-61.3 


60.2-61.6 


tp (year) . . 


2014.6 


2014.5-2014.7 


2014.3-2014.8 



Table 3: Companion dynamical mass and three-dimensional orbit parameters resulting from our MCMC 
analysis. Both 68.2% and 95.4% confidence intervals (CI) are listed. IIR7672B is a benchmark brown dwarf 
with high eccentricity. 



HR7672 is a unique benchmark system because: 
(i) the primary is a GOV star amenable to ultra- 
precise Doppler measurements; (ii) the metallicity 
of the primary is reliably determined from spec- 
troscopy; (iii) the mass and age of the primary are 
well-constrained; and (iv) the parallax has been 
measured with a fractional error of only 1.2%. 

We have monitored IIR7672 with precise RV 
measurements for 24 years, covering 33% of an 
orbital cycle. The observational cadence is high 
compared to other trend stars because this sys- 



tem was also a target of the NASA-UC 77-Earth 



program ([Howard et al.ll2010l) . which began con- 



veniently around the same time as the Doppler 
signal started to show curvature. We have also 
obtained a recent high-contrast image of the sys- 
tem with NIRC2 at Keck, significantly increas- 
ing the astrometric time baseline by extending the 
sky-projected orbital coverage to 13% of a full cy- 
cle. Both data sets have a high («100) signal-to- 
noise ratio and exhibit substantial orbital motion. 
We have performed a joint Doppler-astrometric 
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Fig. 8. — Results from MCMC calculations showing the marginalized posterior distributions of the companion 
period (P), eccentricity (e), inclination (i), argument of periastron (uj), longitude of the ascending node (fi), 
and time of periastron passage (tp). A vertical dashed line indicates the weighted mean value for each 
parameter. Shaded regions correspond to 68.2% and 95.4% confidence intervals respectively. 




Fig. 9. — Companion posterior mass distribution shown before and after taking into account uncertainty in 
the system parallax and mass of the primary star. Uncertainty in the companion mass is limited primarily 
by having only partial orbital phase coverage. 
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Fig. 10. — Probability matrices showing the covariance between the companion mass and aU six orbital 
elements. Each panel shows the posterior distribution of two parameters marginalized over the remaining 
variables including RV offsets and jitter. The contours indicate iso-probability levels corresponding to 68.2% 
(red), 95.4% (light blue), and 99.7% (dark blue) respectively, prior to accounting for uncertainty in the 
system parallax and star mass. A black circle denotes the best-fit values. The companion mass is strongly 
covariant with the orbital period and argument of periastron (and likewise time of periastron passage) . 



MCMC analysis to compute the companion orbital 
elements, true mass, and their uncertainties. 

We find that HR7672B has a near edge-on, i = 
QT.Stoi deg, highly eccentric, e = 0.50t^;J^i, orbit, 
with semimajor axis a — IS.Sj^Q g AU. Presently 
accelerating on its path back towards the primary 
star, both from the perspective of the Earth and 
from the perspective of HR7672A, this companion 
will gradually becoming more difficult to image di- 
rectly. We predict that HR7672B will reach perias- 
tron in the year 2014 just prior to disappearing be- 
hind HR7672A. This companion does not yet have 
JH spectra. We recommend that high-contrast 
imaging programs employing integral-field spec- 
trographs observe this target as soon as possible, 
while the components are reasonably well sepa- 
rated. 

The mass of HR7672B is eS.ltHMj. This 
L4.5 ± 1.5-dwarf therefore resides near the stel- 
lar/substellar boundary, though is likely a brown 



dwarf if not a transition object: 90.1% of the mass 
probability distribution falls below the 72Mj di- 
viding line. Our dynamical measurement is more 
accurate and precise than that obtained previously 
from spectrophotometry and the use of theoreti- 
cal models. With a fractional error of only 4%, it 
is com parable to the results for transiting brow' n 
dwarfs ( Stassun et al ]|2006t I Johnson et alfcoiTal) . 
Uncertainty in the mass will decrease further as 
additional RV and astrometric measurements are 
obtained. 

These observations provide a mass determi- 
nation that is completely independent of the 
companion spectrum, colors, and luminosity (see 
Appendix A for an estimate of the companion 
bolometric luminosity). Presently, few bench- 
mark systems have dynamical mass measure- 
ments and provide such comprehensive infor- 
mation. For instance, GD165B has been well- 
characterized with optical and near-infrared spec- 
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companion mass (Mj^^) 

Fig. 11. — Cumulative summation of the poste- 
rior companion mass distribution. The canonical 
boundary at 72Mj separating stars from brown 
dwarfs is labeled for reference. HR7672B is a 
brown dwarf at 90.1% confidence. Results from 
orbital dynamics are more accurate and precise 
than estimates based on spectrophotometry. 



orbiting a G2V primary. Although RVs have not 
been obtained for the pair, which are nominally 
separated by salOO mas, a near- unity luminosity 
ratio permits a precise individual mass determina - 
tion (|Dupuv et al.ll2009aHKonopackv et al.ll2010l ). 



We conclude that HR7672B is thus a rare and 
precious mass, age, and metallicity benchmark 
brown dwarf companion that may be used to ex- 
plicitly calibrate theoretical evolutionary models 
and synthetic spectral models by anchoring input 
parameters to measured values. 

Detailed substellar companion characterization 
studies are particularly relevant now that the first 
extrasolar planet spectra are securely i n hand 



(jBowler et al.ll2010l: iBarman et all l2011bl lal). As 



the next generation of high- contrast imaging in - 
struments come online (e.g.. iHinklev et al.ll2nilh . 



and more cold bodies orbiting nearby stars are dis- 
covered, HR7672B will serve as a helpful guide for 
the development of more sophisticated theoreti- 
cal models by providing a link between our under- 
standing of stars and planets. 



troscopy, but its age is uncertain because it 



orbits a white d warf ( Kirkpatrick et al. 1999at 



Jones et all I1994I ). HD203030B has also been 



studied across a wide bandpass, conveniently or- 
bits a G8V star, and has a narrow age range 
(130-400 Myr), but its projected separation is 
487 AU, thus prohibiting a dynamical mass 
(|Metchev and Hillenbrandl2006l) . Likewise, HD3651B 
and HD114762B orbit solar- type stars but their 
projected se parations are 480 AU and 128 AU 
respectively (|Liu et al.|[2007t iBowler et al.ll2009l) . 
Many other dynamical measurements have tar- 
geted companions to M-stars , thus helping to cir- 



cumv e nt issues with c ontras t (iZapatero Qsorio et al 
2004 [Ireland et al.l l2008t IDupuv et all l2009bt 



Konopackv et al. 201o t). but M-star metallicities 

are n otoriously difficult to determine ( Johnson and Appj 
20091 ) ■ Low-mass companions to solar- type stars 
have also been studied using Doppler measure- 
ments in combinati on with HST astrometry 
2010l) and Hipparcos intermedi- 



( Benedict et al 



ate astrometry ([Sahlmann et al.l 120111 ). but indi- 
rect observations preclude a comparison between 
spectral models and dynamical mass measure- 
ments by definition. One comparable system, 
HD 130948BC, consists of a double brown dwarf 



6. APPENDIX A: 
OF HR7672B 



THE LUMINOSITY 



We can estimate the luminosity of HR7672B 
using its apparent magnitude and parallax by ap- 
plying a near-infrared bolometric correction de- 
rived from M, L, T field dwarfs. The appar- 
ent magnitudes of HR7672B are measured to be 
TO J « 14.39 ±0.20, TOg = 14.04±0 14, and tuk, = 
13.04±0.10 (|Boccaletti et al.ll2003l) . Of these, the 
K-band magnitude is regarded as most reliable, 
given the contrast ratio between the companion 
and host star. We recalculate the absolute magni- 
tude from lBoccaletti et al.l (|2003l ). owing to the re- 
vised Hipparcos parallax since th e time of publica- 
tion, finding Mk, = 11.79±0.10. iLiu et al.l (|2002t ) 
arrive at a similar answer of Mk^ = 11.8 ± 0.1 us- 
ing relations bet ween spectral type and absolute 
magnitudes from iKirkpatrick et al. I (l2000h . 

To calculate a bolometric magnitude, we next 
derive an estimate for the bolometric correction 



using the empirical relations from lGolimowski et al 
(j2004l ). Correcting for the expected change in 
magnitude between the Kg filter and K filter, 

which results i n an ^ .11 mag increase in 

brightness (i.e., Rudy et al. ( 19961 )). and using a 
spectral-type of L4.5, we find BCk = 3.35 ± 0.13, 
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where the uncertainty is given by the rms scatter 
in the polynomial fit. The bolometric magnitude 
is thus 

Mbol = Mk + BCk = 15.05 ± 0.23. (7) 

Using the Sun as a reference with Mbol,q = 4.74, 
we find the luminosity of HR7672B is equal to 
L = (7.5 ± 1.6) X IO-^Lq. With an age of ^age = 
2.5 ±1.8 Gyr, this result appears t o be consistent 



with theoretical cooling models ([Burrows et al 



19971 ). Applying the Stefan-Boltzmann equation, 
we estimate the radius to be i? « 1.0 ± 0.4i?j as- 
suming an effective temperature in the middle o f 
the 1510-1850 K range found bv lLiu et al.l (|2002l ) . 



7. APPENDIX B: THE SPIN AXIS OF 
HR7672A 

We can also estimate the inclination of the pri- 
mary star rotation axis. Modeling of HR7672A's 
spectrum yields a projected rotational velocity of 
Vsini ~ 2.1 ± 0.7 km/s. Using our measurement 
of logi?^^ = -4.854 ± 0.025 for the chromo- 
spheric Ca I I H a nd K emission, we employ the 
Noves et al.l (jl98J) empirical relations to estimate 
the rotational period, Prot- With B — V = 0.61 
and a convective turnover time of log(r/days) = 
0.99 ± 0.06, we find Prot = 17-5 ± 2.3 days, as- 
suming a mixing length to scale height rat io of 
1.9 (see lWright et al.lbooi iNoves et al.lll984L and 
references therein for details). Combining this re- 
sult with our direct radius measurement (j i2.3.ip . 
we find that the inclination of the stellar rotation 
axis is ij-Qi = 39° ± 16°, where ipot = corre- 
sponds to a pole-on orientation. Thus, the spin 
axis of HR7672A and the orbit of HR7672B are 
not aligned. 
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Dewar Telescope RV oflFset (m/s) RV jitter (m/s) 



6 Lick 24.311:^ 

8 Lick 28.4l^-^ 9.8+2'7 

18 Lick 38.8i^-^ S.Stj'g 

24 Lick ^-'^-I't 

39 Lick 46.9t^;° 9.3;^J 

102 Keck 50. 7t^;^ 

103 Keek G^.Oll'M! 



Table 5: RV offsets and jitter values derived from our MCMC analysis (68.2% confidence interval). 
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Table 6 

DoppLER Radial Velocity Data 



JD-2,440,000 RV (m/s) error (m/s) Instrument Dewar 



7047 71 Q4 


270.71 


0.81 


L 


5 


7373 8210 


232 88 


0.38 


L 


6 


737H Q403 


222.46 


0.68 


L 


g 


7374 8270 


208 96 


0.63 


L 


6 


7374 Qfi31 


241.49 


0.87 


L 


5 


73Q4 7872 


223.48 


1.29 


L 


5 


7430 6837 


260 22 

^ Kjyj * 


1.02 


L 


6 


7431 6884 


273.80 


1.10 


L 


6 


7578 0606 


247.64 


1.29 


L 


6 


7710 8762 


232.07 


1.29 


L 


5 


7846 7068 


201.73 


1.03 


L 


5 


8018 Q8Q5 


191.07 


1.64 


L 


6 


8019 Q631 


216.91 


1.63 


L 


6 


81 1 3 8532 


150.19 


1.02 


L 


5 


8200 61 34 


176.47 


1.03 


L 


5 


8375 00Q3 


158.68 


1.63 


L 


5 


8437 Q603 


173.20 


1.53 


L 


5 


8834 8294 


133.43 


0.75 


L 


8 


8847 8253 


128.49 


0.93 


L 


8 


8905 7100 


136.85 


0.84 


L 


8 


9123 0038 


91.09 


1.38 


L 


8 


9171 8390 


103.59 


1.11 


L 


8 


91 73 8394 


98.98 


1.26 


L 


8 


91 74 9026 


1 1 75 


1.12 


L 


8 


9200 8275 


100.12 


1.04 


L 


8 


Q970 7Q7Q 


1 0^ nn 


1 OS 


T 


Q 
O 


9280.7483 


97.92 


0.99 


L 


8 


9587.7636 


80.56 


0.98 


L 


8 


9589.7614 


82.45 


1.03 


L 


8 


9602.7461 


101.85 


0.96 


L 


8 


9622.7247 


90.93 


1.26 


L 


8 


9680.5946 


66.84 


0.76 


L 


39 


9872.9322 


54.42 


0.49 


L 


39 


9893.8756 


61.43 


0.50 


L 


39 


9913.9177 


45.84 


0.51 


L 


39 


9984.6846 


56.55 


0.51 


L 


39 


10215.8672 


8.85 


0.78 


L 


39 


10263.8492 


26.16 


0.54 


L 


39 


10299.8286 


33.72 


0.61 


L 


39 


10304.8468 


18.14 


0.56 


L 


39 


10601.8849 


0.00 


0.01 


K 


102 
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Table 6 — Continued 



JD-2,440,000 RV (m/s) error (m/s) Instrument Dewar 



10602.9911 


3.78 


1.29 


K 


102 


10604.0217 


3.69 


1.19 


K 


102 


10605.0298 


2.92 


1.30 


K 


102 


10606.0621 


-5.73 


1.30 


K 


102 


10606.9658 


0.77 


1.38 


K 


102 


10607.9157 


3.79 


1.37 


K 


102 


10613.9009 


-2.26 


0.79 


L 


39 


10614.9313 


-0.37 


0.75 


L 


39 


10640.8998 


6.26 


0.75 


L 


39 


10656.8269 


7.41 


0.76 


L 


39 


10793.5628 


1.76 


0.51 


L 


39 


10979.9278 


-6.21 


1.16 


L 


18 


11048.8474 


5.92 


1.17 


L 


18 


11062.7690 


-9.13 


1.23 


L 


18 


11068.8450 


-23.71 


1.45 


K 


102 


11069.8909 


-29.20 


1.34 


K 


102 


11070.9108 


-19.72 


1.31 


K 


102 


11071.8471 


-13.15 


1.60 


K 


102 


11072.8352 


-15.16 


1.34 


K 


102 


11074.8549 


-22.96 


1.31 


K 


102 


11075.7848 


-22.30 


1.18 


K 


102 


11300.0046 


-33.20 


0.78 


L 


39 


11303.9778 


-33.21 


0.80 


L 


18 


11305.9730 


-32.72 


0.88 


L 


18 


11311.1002 


-41.88 


1.23 


K 


102 


11312.0886 


-40.53 


1.37 


K 


102 


11313.1019 


-33.60 


1.39 


K 


102 


11314.1248 


-39.07 


1.04 


K 


102 


11367.9059 


-53.11 


1.51 


K 


102 


11368.8710 


-49.16 


1.72 


K 


102 


11370.0152 


-45.70 


1.28 


K 


102 


11372.0268 


-41.96 


1.64 


K 


102 


11373.0623 


-46.15 


1.59 


K 


102 


11373.8199 


-51.10 


1.48 


K 


102 


11389.8655 


-51.27 


0.80 


L 


18 


11392.8306 


-42.50 


0.89 


L 


18 


11410.8762 


-52.99 


1.52 


K 


102 


11411.8717 


-55.86 


1.39 


K 


102 


11416.7036 


-29.17 


0.88 


L 


18 


11439.8221 


-53.96 


1.46 


K 


102 


11467.6397 


-28.58 


1.57 


L 


18 
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Table 6 — Continued 



JD-2,440,000 RV (m/s) error (m/s) Instrument Dewar 



11468.6072 


-34.72 


1.31 


L 


18 


11679.0915 


-39.86 


1.85 


K 


102 


11703.0748 


-46.32 


1.50 


K 


102 


11751.8211 


-50.03 


0.88 


L 


18 


11754.9418 


-49.21 


1.33 


K 


102 


11793.8250 


-53.33 


1.63 


K 


102 


11815.6476 


-45.26 


1.48 


L 


18 


12031.0398 


-79.77 


1.53 


K 


102 


12041.9946 


-78.96 


1.22 


L 


18 


12075.9299 


-81.23 


0.98 


L 


18 


12099.0483 


-87.46 


1.60 


K 


102 


12115.8858 


-87.66 


0.95 


L 


18 


12117.8641 


-82.71 


0.82 


L 


18 


12127.9831 


-100.94 


1.83 


K 


102 


12133.8002 


-98.58 


1.62 


K 


102 


12140.7744 


-69.50 


0.87 


L 


18 


12181.6483 


-84.26 


0.99 


L 


18 


12186.6687 


-75.66 


1.15 


L 


18 


12188.7009 


-77.17 


1.29 


L 


18 


12202.6766 


-91.70 


1.09 


L 


18 


12391.1392 


-120.49 


1.91 


K 


102 


12508.7589 


-120.32 


1.18 


L 


24 


12509.7462 


-134.04 


1.14 


L 


24 


12515.7825 


-115.41 


1.12 


K 


102 


12575.6912 


-116.79 


1.69 


K 


102 


12848.9031 


-131.74 


1.78 


K 


102 


12855.9988 


-136.42 


1.67 


K 


102 


12894.7297 


-139.06 


1.45 


L 


24 


13181.0272 


-171.06 


1.48 


K 


102 


13203.8281 


-160.71 


1.90 


L 


24 


13303.7113 


-172.17 


0.84 


K 


103 


13544.9448 


-203.04 


1.47 


L 


24 


13551.0081 


-195.05 


0.55 


K 


103 


13556.9127 


-190.72 


1.37 


L 


24 


13561.8937 


-201.13 


1.58 


L 


24 


13589.8255 


-186.90 


2.81 


L 


24 


13603.9205 


-195.75 


0.74 


K 


103 


13926.0386 


-212.10 


0.67 


K 


103 


13926.8957 


-217.61 


1.65 


L 


24 


13954.7593 


-226.25 


1.97 


L 


24 


13958.7828 


-221.32 


1.29 


L 


24 



24 



Table 6 — Continued 



JD-2,440,000 


RV (m/s) 


error (m/s) 


Instrument 


Dewar 


13982.8138 


-220.35 


0.62 


K 


103 


13989.8700 


-224.79 


3.04 


L 


24 


14337.0434 


-235.41 


1.44 


K 


103 


14546.1518 


-240.05 


1.22 


K 


103 


14671.9587 


-240.90 


1.21 


K 


103 


14673.9656 


-248.17 


1.16 


K 


103 


14675.8841 


-247.46 


0.93 


K 


103 


14688.8742 


-239.27 


1.47 


K 


103 


14689.9519 


-241.90 


1.36 


K 


103 


14697.9197 


-222.74 


2.94 


L 


24 


14717.9307 


-244.81 


1.25 


K 


103 


14718.9864 


-246.41 


1.27 


K 


103 


14719.9911 


-239.01 


1.18 


K 


103 


14720.9566 


-225.02 


1.22 


K 


103 


14721.9701 


-218.99 


1.19 


K 


103 


14722.8696 


-226.22 


1.20 


K 


103 


14724.9263 


-233.32 


1.27 


K 


103 


14725.8416 


-232.83 


1.23 


K 


103 


14726.9581 


-231.17 


1.16 


K 


103 


14727.8365 


-236.50 


1.21 


K 


103 


14777.8159 


-235.80 


1.43 


K 


103 


14808.7068 


-234.60 


1.35 


K 


103 


14929.1287 


-234.40 


1.43 


K 


103 


14930.1221 


-237.22 


1.34 


K 


103 


14935.1295 


-230.98 


1.22 


K 


103 


1 SflzLI S7ni 
iOU^l .O / Ul 


99Q <^<^ 
-zzy .00 


1 '?zL 
± .o^ 


K 


lUO 


15076.7430 


-230.53 


1.29 


K 


103 


15079.7372 


-232.51 


1.34 


K 


103 


15106.9076 


-219.21 


1.36 


K 


103 


15169.7567 


-230.30 


0.88 


K 


103 


15290.1483 


-218.21 


1.36 


K 


103 


15313.1369 


-224.11 


1.32 


K 


103 


15319.0997 


-226.52 


1.48 


K 


103 


15404.9095 


-214.63 


1.35 


K 


103 


15437.0096 


-206.43 


1.35 


K 


103 


15542.6875 


-198.94 


1.15 


K 


103 


15723.0122 


-163.58 


1.38 


K 


103 


15785.9666 


-169.26 


1.58 


K 


103 


15842.7963 


-158.36 


1.57 


K 


103 
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Note. — Raw RV measurements used for orbital analysis. Uncer- 
tainties correspond to the error due to photon noise only. Observa- 
tions obtained at Lick and Keck are labeled in the instrument column. 
Several different hardware configurations have been used over the past 
24 years. They arc labeled according to dewar number. Each require 
different RV offset and jitter values. 
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Table 7 

Literature Photometry for HR7672 



Filter 




¥ Oj± L_LV_. \ XXXClJtL 1 


Reference 


Johnson 


V 


5.80 ± 0.05 


Johnson et al. 


(1966) 


Johnson 


B 


6.410 ± 0.05 


Johnson et al. 


(1966) 


Johnson 


U 


6.500 ± 0.05 


Johnson et al. ('19661 


Johnson 


V 


5.80 ± 0.05 


Cowlev et al. 


(1967) 


Johnson 


B 


6.390 ± 0.05 


Cowlev et al. 


(1967) 


Johnson 


U 


6.500 ± 0.05 


Cowlev et al. 


(1967) 


Johnson 


V 


5.80 ± 0.05 


Johnson and Knuckles fl957) 


Johnson 


B 


6.410 ± 0.05 


Johnson and Knuckles fl957') 


Johnson 


U 


6.500 ± 0.05 


Johnson and Knuckles ('19571 


Johnson 


V 


5.77 ± 0.05 


Niconov et al. 


(1957) 


Johnson 


B 


6.360 ± 0.05 


Niconov et al. 


(1957) 


Johnson 


V 


5.79 ± 0.05 


MermiUiod ( 1986) 


Johnson 


B 


6.390 ± 0.05 


Mermilliod (1986) 


Johnson 


U 


6.490 ± 0.05 


Mermilliod ( 


1986) 


Johnson 


V 


5.799 ± 0.05 


Moffett and Barnes { 1980) 


Johnson 


B 


6.402 ± 0.05 


Moffett and Barnes { 1980) 


Johnson 


V 


5.80 ± 0.05 


Nicolet (1978) 


Johnson 


B 


6.410 ± 0.05 


Nicolet C1978) 


Johnson 


U 


6.500 ± 0.05 


Nicolet (1978) 


Johnson 


U 


6.5 ± 0.05 


Barnes et al. 


(1997) 


Johnson 


u 


6.5 ± 0.05 


Barnes et al. 


(1997) 


Johnson 


u 


6.53 ± 0.05 


Barnes et al. 


(19971 


Johnson 


u 


6.5 ± 0.05 


Barnes et al. 


(1997) 


Johnson 


u 


6.49 ± 0.05 


Barnes et al. 


(1997) 


Johnson 


B 


6.41 ± 0.05 


Barnes et al. 


(1997) 


Johnson 


B 


6.39 ± 0.05 


Barnes et al. 


(1997) 


Johnson 


B 


6.41 ± 0.05 


Barnes et al. 


(1997) 


Johnson 


B 


6.41 ± 0.05 


Barnes et al. 


(1997) 


Johnson 


B 


6.36 ± 0.05 


Barnes et al. 


(1997) 


Johnson 


B 


6.394 ± 0.05 


Barnes et al. 


(1997) 


Johnson 


B 


6.39 ± 0.05 


Barnes et al. 


(1997) 


Johnson 


B 


6.402 ± 0.05 


Barnes et al. 


(1997) 


Johnson 


V 


5.8 ± 0.05 


Barnes et al. 


(19971 


Johnson 


V 


5.8 ± 0.05 


Barnes et al. 


(19971 


Johnson 


V 


5.79 ± 0.05 


Barnes et al. 


(19971 


Johnson 


V 


5.8 ± 0.05 


Barnes et al. 


(19971 


Johnson 


V 


5.77 ± 0.05 


Barnes et al. 


(19971 


Johnson 


V 


5.797 ± 0.05 


Barnes et al. 


(1997) 


Johnson 


V 


5.79 ± 0.05 


Barnes et al. 


(1997) 


Johnson 


V 


5.799 ± 0.05 


Barnes et al. 


(1997) 


Stromgren 


u 


7.681 ± 0.08 


Hauck and Mermilliod (1998) 
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Table 7 — Continued 



Filter 


Value (mag) 


Reference 


Stromgren v 


6.765 ± 0.08 


Hauck and Mermilliod f 19981 


Stromgren b 


6.187 ± 0.08 


Hauck and Mermilliod fl998) 


Stromgren y 


5.80 ± 0.08 


Hauck and Mermilliod fl998~) 


Stromgren u 


7.673 ± 0.08 


Olsen (1994) 


Stromgren v 


6.757 ± 0.08 


Olsen (1994) 


Stromgren b 


6.184 ± 0.08 


Olsen (1994) 


Stromgren y 


5.80 ± 0.08 


Olsen (1994) 


Stromgren u 


7.675 ± 0.08 


Fabregat and Reglero (1990) 


Stromgren v 


6.761 ± 0.08 


Fabrcgat and Reglero (1990) 


Stromgren b 


6.184 ± 0.08 


Fabrcgat and Reglero fl990") 


Stromgren y 


5.80 ± 0.08 


Fabregat and Reglero C1990~) 


Stromgren u 


7.682 ± 0.08 


Olsen 


(1993) 


Stromgren v 


6.761 ± 0.08 


Olsen 


(1993) 


Stromgren b 


6.187 ± 0.08 


Olsen (1993) 


Stromgren y 


5.80 ± 0.08 


Olsen 


a993) 


Stromgren u 


7.682 ± 0.08 


Crawford et al. fl966) 


Stromgren v 


6.775 ± 0.08 


Crawford et al. (1966") 


Stromgren b 


6.189 ± 0.08 


Crawford et al. (1966) 


Stromgren y 


5.80 ± 0.08 


Crawford et al. (1966) 


Stromgren u 


7.682 ± 0.05 


Olsen 


(1994) 


Stromgren u 


7.695 ± 0.05 


Olsen 


C1994) 


Stromgren u 


7.681 ± 0.05 


Olsen 


C1994) 


Stromgren u 


7.672 ± 0.05 


Olsen 


C1994) 


Stromgren b 


6.189 ± 0.05 


Olsen 


C1994) 


Stromgren b 


6.204 ± 0.05 


Olsen 


C1994) 


Stromgren 6 


6.186 ± 0.05 


Olsen 


(1994) 


Stromgren b 


6.183 ± 0.05 


Olsen 


(1994) 


Stromgren v 


6.775 ± 0.05 


Olsen 


C1994) 


Stromgren v 


6.781 ± 0.05 


Olsen 


C1994) 


Stromgren v 


6.76 ± 0.05 


Olsen 


(1994) 


Stromgren v 


6.756 ± 0.05 


Olsen 


C1994) 


Stromgren y 


5.8 ± 0.05 


Olsen 


C1994) 


Stromgren y 


5.82 ± 0.05 


Olsen 


C1994) 


Stromgren y 


5.799 ± 0.05 


Olsen 


(1994) 


Stromgren y 


5.799 ± 0.05 


Olsen 


(1994) 


Geneva V 


5.794 ± 0.08 


Rufener (19761 


Geneva VI 


6.541 ± 0.08 


Rufener (1976) 


Geneva B 


5.571 ± 0.08 


Rufener (1976) 


Geneva Bl 


6.651 ± 0.08 


Rufener (1976) 


Geneva B2 


6.902 ± 0.08 


Rufener (1976) 


Geneva U 


6.888 ± 0.08 


Rufener (1976) 
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Table 7 — Continued 



Filter 


Value (mag) 


Reference 




Geneva G 


6.846 ± 0.08 


Rufener (1976) 




2Mass Ks 


4.388 ± 0.0230 


Cutri et al. C20031 




DDO 48 


6.017 ± 0.05 


McClure and Forrester 


fl98D 


DDO 45 


7.02 ± 0.05 


McClure and Forrester 


fl981) 


DDO 42 


7.643 ± 0.05 


McClure and Forrester 


fl981^ 


DDO 41 


7.658 ± 0.05 


McClure and Forrester 


flQSl) 


DDO 38 


6.863 ± 0.05 


McClure and Forrester 


fl981) 


DDO 35 


7.767 ± 0.05 


McClure and Forrester 


fl981) 
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